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Abstract 


We use five and one-haif years of limb- and nadir-viewing temperature mapping obser- 
vations by the Composite Infrared Radiometer-Spectrometer (CIRS) on the Cassini 
Saturn orbiter, taken between July 2004 and December 2009 (L$ from 293° to 4° ; 
northern mid-winter to just after northern spring equinox), to monitor temperature 
changes in the upper stratosphere and lower mesosphere of Titan. The largest changes 
are in the northern (winter) polar stra.topa.use, which has declined in temperature by 
over 20 1< between 2005 arid 2009. Throughout the rest of the mid to upper strato- 
sphere and lower mesosphere, temperat ure changes are less than 5 K. In the southern 
hemisphere, temperatures in the middle stratosphere near 1 mbar increased by 1 to 
2K from 2004 through early 2007, then declined by 2 to 4K throughout 2008 and 
2009, with the changes being larger at more polar latitudes. Middle stratospheric 
temperatures at mid-northern latitudes show a small 1 to 2 1\ increase from 2005 
through 2009. At north polar latitudes within the polar vortex, temperatures in the 
middle stratosphere show a ~ 4 1< increase during 2007, followed by a comparable 
decrease in temperatures in 2008 and into early 2009. The observed temperature 
changes in the north polar region are consistent with a weakening of the subsidence 
within the descending branch of the middle atmosphere meridional circulation. 


Key words: Titan, atmosphere. Atmospheres, structure. Atmospheres, dynamics, 
Inf r ared O bs ervati ons 
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i 1. Introduction 


2 With an obliquity of approximately 26.7° , there are strong seasonal variations of 

3 the insolation of Titan’s atmosphere. Although the radiative timescales in the tropo- 

4 sphere are longer than the Titan year (Smith et ah, 1981), in the mid-stratosphere 
s and above the radiative response time becomes short enough that seasonal variations 
e in temperature, and thus dynamics, are expected (Fiasar et ah, 1981). 

r The earliest measurements of Titan’s atmospheric temperatures were from the 
a Voyager 1 infrared spectrometer IRIS, which observed Titan in early northern spring 

9 ( Ls ~ 9°). Retrievals of temperatures at 0.4 and 1 mba.r by Fiasar and Conr&th (1990) 

10 showed a hemispheric asymmetry, with temperatures near 60° N approximately 15 K 
n colder than equatorial and low southern latitudes, with a smaller drop of about 5 K 
is at 60° S. This temperature asymmetry is somewhat larger than expected given ra- 
i 3 diative cooling timescales at these levels of about 1 Earth year (one Titan year is 
u 29.-5 Earth years). Fiasar and Conrath (1990) proposed that, because of the coupling 
is of the temperatures and winds through gradient wind balance, the temperatures re- 
16 spond to the solar forcing on the longer dynamical timescale on which the meridional 
I? circulation transports zonal momentum between hemispheres. Alternately, Bezard 
is et al. (1995) suggested that the temperature asymmetry is the result of hemispheric 

1 9 asymmetries in the opacity sources responsible for solar heating and radiative cooling. 

20 Stratospheric temperature retrievals from Cassini Composite Infrared Spectrometer 

2 1 (CIRS) data taken around northern mid-winter (Fiasar et ah, 200-5; Achterberg et ah, 

22 2008a), showed a somewhat stronger asymmetry, with 1 mhar temperatures poleward 

23 of 60° N 20-30 K colder than the equator. The most remarkable result from the CIRS 
a* temperature data is the discovery that the stratopause was over 20 K warmer and 
25 over a scale height in altitude higher at high northern winter latitudes than at the 
so equator (Achterberg et ah, 2008a). As the winter stratopause is high enough in al- 
27 titude to be above the polar shadow, and there is an increased haze opacity in the 
2 * winter polar hood (Smith et al, 1981), increased radiative heating may contribute to 
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the warm stratopause. Another possibility suggested by Achterberg et al. (2008a) is 
30 adiabatic heating in the descending branch of the meridional circulation. Subsidence 
3 ! has also been invoked to explain the enhanced abundance of nitriles and some hydro- 

32 carbons seen at high northern latitudes by Voyager and Cassini (Rannou et ah, 2004; 

33 Lebonnois et ah, 2009; Teanby et ah, 2009). 

34 With five and a half years of data now available from Cassini, covering just over 

35 two- thirds of a season from early northern winter to just after northern spring equinox 
3 f> (L$ from 293° to 4°), we have used both nadir and limb viewing temperature map- 

37 ping observations from CIRS to look for temporal changes in temperature in the tapper 

38 stratosphere and lower mesosphere. Section 2 describes the data sets used and the 
m temperature retrieval procedure. Section 3 presents the observed time variations of 

40 the retrieved temperatures and of the zonal mean gradient winds consistent with the 

4 1 temperatures. Section 4 compares the observations with the earlier Voyager observa- 

42 tions and with general circulation models of Titan, and discusses the implications of 

43 the observations for the mean meridional circulation. 

44 2. Observations and Data Analysis 

45 To obtain temperatures in the upper stratosphere, we use thermal infrared spec- 

46 tra from focal plane 4 (FP4) of CIRS. which covers the spectral range from 1100 

47 to 1400 cm -1 with a spectral resolution adjustable between approximately 0.5 and 
4 s 15.5 cm -1 . FP4 is a 10-element linear array, each pixel of which has a square field, of 
49 view of 0.28 mrad. The instrument and its operation is described in detail by Flasar 
so et al. (2004). The data used cover the time period from the TO flyby on 2 July 2004 
si (L$ — 293°). just after orbit insertion, through the T63 flyby on 11 December 2009 

52 (Ls = 4°), just after Titan’s northern spring equinox. 

5 3 Two types of temperature mapping observations are used: nadir-viewing maps 

54 which usually cover the visible hemisphere of Titan, and limb-viewing maps which 

55 provide better altitude coverage than nadir-viewing observations, but which cover 
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56 only a single longitude for each observation. For the nadir-viewing maps, a series of 
m continuous slews, perpendicular to the 1x10 FP4 array, are used to cover the visible 
5s disk from a spacecraft distance of 250.000 to 400,000 km, giving a spatial resolution 
59 of 1.5 s to 2.5 s of arc at, disk center, and the data is taken at a spectral resolution of 
so 2.8 cm" 1 . A list of the nadir observations is given in Table 1. For the limb- viewing 
si maps, the array is placed perpendicular to the limb of Titan, with two separate 
ea vertical placements of the array, one above the other with an approximately 3 pixel 
ss overlap. Data, are taken every 5° in latitude over a roughly 90° latitude range. At 
«4 each position of the array, 10 to 12 spectra are acquired in each detector at a spectral 
c-5 resolution of 15.5 cm" 1 , and are averaged together to improve the signal- to-noise ratio, 
ee Spectra taken during the slews between array positions are ignored. The limb maps 
5- are taken at a range of about 120,000 km, giving a vertical spatial resolution of about, 
as 40 km, somewhat less than the pressure scale height in the stratosphere. A list of the 
63 limb observations used is given in Table 2. The latitude and time coverage provided 
M by both sets of observations is shown in Fig. 1 . 

n Temperatures are retrieved from the observed spectra using the procedure de- 

72 scribed in detail by Achterherg et al. (2008a); only a very brief summary is provided 

73 here. For the nadir maps, temperatures are retrieved from individual spectra by 

74 a constrained linear inversion algorithm (Conrath et al., 1998). using the spectral 
n range from 1251 to 1311cm" 1 within the P and Q branches of the v 4 band of CH 4 , 
76 The forward model uses 200 layers, equally spaced in log-pressure between 1456 and 
?7 0.001 mbar. CH 4 iq band absorption is calculated using the correlated-^ approxima- 
te tion (Lacks and Oinas, 1991), with line data taken from the 2003 release of the GEISA 
79 line atlas (Jacqninet- Husson et al, 2005), assuming a stratospheric mole fraction of 
so 0.0141 (Niemann et al, 2005). The retrieved temperatures are valid between roughly 
n 0.2 and 8 mbar, except at north polar latitudes where the region of validity shifts 
82 to approximately 0.1 mbar to 2 mbar because of the warmer upper stratosphere and 
as colder middle and lower stratosphere. Outside of the pressure range that is strongly 
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constrained by the data, the retrieved temperatures relax to the initial guess profile. 

For the limb maps, the large vertical extent of the atmosphere relative to the radius 
of Titan requires the use of a 2-dimensional retrieval algorithm (Achterberg et a!., 
2008a), in which all of the data from a single map is used simultaneously to retrieve 
the temperature and aerosol abundance as a function of both latitude and altitude, 
assuming no variations with longitude, as well as altitude pointing corrections. The 
spectral range used 1210 to 1315 cm -1 , The forward model uses 200 layers equally 
spaced in altitude between the surface and 550 km, and latitude points spaced every 
5° , with the latitude boundaries dependent upon the latitude range of the data used. 
Again, CH.j absorption is included using the correlated-^; approximation, and haze 
absorption is also included using a gray absorber with no scattering. The retrieved 
temperatures are valid between approximately 0.005 and 5 mbar, except northwards 
of about 60° N where they are valid between about 0.005 and 2 mbar. For both sets 
of retrievals, the initial guess for the retrievals was the latitude-dependent, smoothed, 
time- and zona liy-aver aged temperature cross section of Achterberg et al, (2008a, 
Fig. 9). 

3. Results 

Figure 2 illustrates the large-scale changes in Titan's stratospheric thermal struc- 
ture between the early part of the Cassini mission (northern mid-winter) and the 
period just before northern spring equinox. Cross-sections of zonal mean tempera- 
tures averaged over flybys TO through T13 ( July 2004 through March 2006) are shown 
in the top panel, and over T45 through T63 (July 2008 through December 2009) in 
the bottom panel. The averages were constructed by averaging both the nadir- and 
limb-viewing data in 5° latitude bins, and smoothing the resulting averages by 3 ap- 
plications of a. sliding 10° rectangular window. The middle panel shows the change in 
temperature between the early and late temperature averages. The largest changes 
in temperature occur at around 60° N, near the northern boundary of the winter po- 
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111 lar vortex, with only small changes in the southern hemisphere. The most prominent 
m change is a strong cooling of the warm winter polar stratopause by as 20K, along with 
1,3 a general cooling of the stratopause and lower mesosphere. In the mid-stratosphere, 
m temperature changes are smaller than in the lower mesosphere, with warming of 2-3K 
ns in the northern hemisphere and weak cooling in the southern hemisphere. 

ns 3.1. Limb Maps: Lower Mesosphere and Stratopause 

i,7 Figure 3 shows a time sequence of temperature cross-sections from individual limb 
ne maps with northern hemisphere coverage. The dashed line on the figure indicates the 
H9 stratopause, calculated as the pressure where the temperature is a local maximum at 

120 each latitude. The warm, elevated north polar stratopause has been cooling from 2005 

121 through 2009 at a rate of about 5 K per year. In addition, the transition between 

122 the ~ 0.1 mbar near-equatorial stratopause and the elevated polar stratopause has 
m become much more abrupt, and in the T59 and T63 maps (Ls — 359° and L,$ = 
-.24 4°). the stratopause has essentially disappeared at northern midlatitudes, with the 
os temperatures profile becoming nearly isothermal for pressures less than about 0.3 
ns mbar. 

127 Figure 4 shows a time sequence of equatorial temperature cross-sections from 
we individual limb maps. From 2005 through 2007, changes to the temperature structure 

129 are small. From late 2007 through mid 2009, temperatures in the lower mesosphere 

130 between about 0.1 and 0.01 mbar decrease by approximately 5 K, accompanied by a 
ns shift in the stratopause to higher pressures. 

is: Figure 5 shows a time sequence of southern hemisphere temperature cross-sections 

ns from individual limb maps. Coverage of high southern latitudes from the limb maps is 

134 unfortunately sparse (see Fig. 1). Little evidence of temperature changes can be seen, 

135 apart from indications of cooling in the lower mesosphere at higher latitudes, and the 
os disappearance of the slight elevation of the stratopause at south polar latitudes. 
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ij7 3.2. Nadir Maps: Middle Stratosphere 

ns To examine the temporal variations in the middle stratosphere, temperatures re- 
us trieved from the nadir maps for each flyby were zonally averaged in 5° wide latitude 
bins. In calculating the averages, the variation of latitude along the observed ray 
mi paths, which is significant because of the large vertical extent of Titan’s atmosphere, 
542 was taken into account. Figure 6 shows a time series of the zonal mean tempera, - 
i4j tures at 1 mbar for a representative set. of latitudes. Unfortunately, flyby geometries 
144 did not allow for good nadir views of latitudes poleward of approximately 60° N or 
us 70° S for the first two years of the mission, except for the south pole on the TO flyby 
i4o immediately after orbit, insertion. 

547 The north polar data at; 75° N and 85° N indicate a temperature increase of ~ 3 K 
Ha between early 2007 and early 2008, but a similar decrease from early 2008 through 

149 spring equinox, despite increasing sunlight in the north polar region as spring ap- 

150 proaches. The mid- northern hemisphere temperatures (30° N - 60° N) show consider- 
5bi able scatter, but indicate a temperature increase of 2 to 3K from 2005 to 2009, con- 
152 sistent with the larger time averages shown in Fig. 2. Near equatorial temperatures 
553 are constant over the time of the observations. At mid-to-high southern latitudes, 

154 temperatures increased by about 1 K between mid- 2004 and early 2007; there is then 

155 clear temperature decrease beginning in late 2007, with the amplitude of the decrease 

156 from late 2007 to mid- 2009 varying from ~ 1 K around 45° S to ~ 4.K near the south 
is? pole. 


158 3.3. Zonal Winds 


is 9 Using the assumption of gradient wind balance - that the meridional pressure 

i eo gradient is balanced by the horizontal component of the the sum of the Coriolis and 

161 centrifugal forces, which is expected to hold for Titan - the zonal wind velocity u is 

162 related to the meridional gradient of the temperature T by 
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163 Here 0 is Titan’s rotation frequency, 0 is latitude, g is gravitational acceleration, 
and for a thick atmosphere, the “vertical’’ derivative dz$ is taken along the direction 
i65 parallel to the rotation axis (Fiasar et al., 2005). To minimize the noise from taking 
iw derivatives of the temperature field, and to get full latitude coverage. Eq. 1 is applied 
i57 to the averaged temperature fields of Fig. 2. Following Achterberg et al. (2008a), we 

168 use as a boundary condition a wind in solid body rotation at 10 mhar with an angular 
>69 velocity of four times the solid body rotation rate, consistent with the winds at that 
no level measured by the Huygens Doppler Wind Experiment (Bird et al., 2005; Folkner 
Hi et ah, 2006) . 

ira The resulting winds are shown in Fig. 7. along with the change in the wind between 
1-3 the early and late averages. Changes in the southern hemisphere winds are small, 
174 as expected from the small changes in the southern hemisphere temperatures. The 
its main change in the winds between 2005 and 2009 is an extension of the strong winter 
ns hemisphere jet to higher altitudes; the weakening meridional temperature gradients at 

177 the latitude of the jet correspond to a slower decay of the jet with increasing altitude. 

178 $.4. Stratospheric Pole Offset 

179 Analysis of mid-stratospheric temperatures by Achterberg et al. (2008b) showed 
iso that, zonal variations of temperature were dominated by a zonal wavenumber 1 feature 
lai which was stationary in either a solar-fixed or stellar-fixed reference frame to within 

183 the observational uncertainty. The meridional structure of the wavenumber 1 feature 
is3 was consistent with Titan’s stratospheric temperatures and winds being symmetric 

184 about an axis that is offset from the then-current IAU definition of Titan’s solid body 

185 rotation axis by 4.1° ± 0.2°, with the north pole tilted towards a direction 76° dfc T 
m west of the subsolar longitude. Subsequently, similar offsets were discovered in the 
is7 symmetry axis of the hemispheric north-south haze asymmetry (Roman et al., 2009) 
lea and of the HCN abundance (Teanby et al., 2010). 

169 Recent GCM simulations by Tokano (2010) produce a stratospheric atmospheric 
angular momentum (AAM) vector which is offset from the rotation pole by a few 
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j m degrees and which precesses westward with a period of one Titan solar day, as a 

| 192 result of forcing by solar thermal tides. These results are broadly consistent with the 

i» observations, except that the AAM vector in the GCM simulations is tilted toward the 
194 anti-solar longitude in northern winter. Tokano (2010) also found that the angle of the 
sat AAM tilt in the GCM varies strongly with season, with a maximum tilt amplitude 
i66 occurring after the solstices and the tilt becoming near zero shortly after equinox, 
m with the tilt angle in the mid-stratosphere dropping by ~30% to 50% over the period 
las of the G’IRS observations. 

199 To look for temporal variations of the stratospheric tilt, the procedure described 

200 by Achterberg et al. (2008b) for fitting the symmetry axis of the temperature field was 

201 applied to the retrieved 1 mbar temperature maps for each Titan flyby. The results 

202 are shown in Fig. 8 in both sun-fixed and inertial (star-fixed) reference frames. The 
208 error bars shown are 1 -a uncertainties resulting from uncertainties in the retrieved 
w temperatures. Navigation of support imaging for the earliest Titan temperature maps 
20$ indicated uncertainties in the reconstructed pointing data of one-third to one-half of 
a* a CIRS FP4 field of view, and retrieved pointing corrections from the limb retrievals 

207 are also in general one-third of a field of view or less, corresponding to errors of about 

208 1° of arc on Titan. As the pointing errors will be moderately strongly correlated over 
303 a temperature mapping observation, it is likely that much of the considerable scatter 
2:o in Fig. 8 is from pointing errors. 

ns Despite the scatter, the amplitude of the fitted pole offset shows no indication of 

212 having decreased over the period of the observations as predicted by the model of 

213 Tokano (2010); if anything the data suggest a slight increase in the amplitude of the 
aw tilt. Linear fits to the azimuth angle of the tilt as a function of time, shown as dotted 
sis lines in Fig. 8, suggest that the tilt is closer to stationary in the inertial reference 
2i6 frame than in the solar reference frame; the fitted long-term drift rates of the azimuth 
as? are (29 ± 8) x 10 -S degs _1 in the sun-fixed frame and (3.5 ± 7.9) x 10~ 8 deg s' 1 in 
2i$ the inertial frame. The plots of tilt azimuth also show an apparent oscillation with 
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219 a period of about 3 years. Calculation of a Lomb-Seargle periodogram (Press and 

220 Teukolsky, 1988) from the tilt azimuth data in the inertial frame gives a maximum 
25 i signal at a period of 2.73 Earth years, with a false alarm probability of 11.5% (i.e. 
H 2 statistically significant at the 88.5% level). With the large scatter in the tilt data, 

223 and the correspondingly low statistical significance of the results, further data to 

224 be obtained during the Cassini Solstice Mission are needed to clearly determine the 
2 » behavior of the stratospheric pole offset. 

226 4. Discussion 

227 JJ- Comparison with Voyager 1 IRIS 

258 Prior to Cassini, the only measurements of meridional temperature profiles on 
229 Titan were from the Infrared Interferometer Spectrometer (IRIS) on Voyager 1 (Flasar 
2 M et al., 1981; Flasar and Conrath. 1990; Coustenis and Bezard, 1995), which flew by 

23 1 Titan in November 1980. during early northern spring. (L$ = 9°). Figure 9 show's the 

232 0.4 mbar and 1.0 mbar temperatures retrieved from Voyager 1 IRIS by Flasar and 

233 Conrath (1990), along with temperatures at the same pressures retrieved from CIRS 

234 data from flybys T45 through T59 (July 2008 to July 2009, Lg =■ 347° - 359°), taken 
ass slightly less than one Titan year after the IRIS data. Overall, the two data sets are 

236 consistent with each other, with two main differences. First, the Voyager data show 

237 an indication of the meridional temperature gradient becoming steeper at 60° S, which 

238 is not seen in the CIRS data, but which is consistent with the the Voyager data being 
aw taken in early southern spring, while the Cassini data is from late southern winter. 
243 Secondly, the Voyager data are somewhat colder at the equator at 0.4 mbar than the 
24s CIRS data. 

242 4.2. Comparison with Titan General Circulation Models 

243 The primary objective of Titan general circulation modeling is to contribute to an 

244 understanding of the processes that control the structure and dynamics of the atmo- 
2 « sphere. This modeling objective can only be met through the comparison of model 
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predictions with the available measurements of the spatial structure and temporal 
evolution of atmospheric temperatures, wind fields, gas abundances, and aerosols, 
3« By adjusting parameterizations and systematically introducing new mechanisms, the 
m relative importance of various processes can be established. The Cassini mission has 

250 obtained a wealth of new information on Titan's atmosphere that can establish im- 

251 portant constraints. Although the fraction of the Titan seasonal cycle covered thus 

252 far by CIRS observations is limited (L$ from 293° through 4°), temporal changes in 

253 the temperature field (Fig. 2) and the associated gradient thermal wind field (Fig. 7) 

254 are clearly evident. These results can now be compared with circulation model pre- 

255 dictions. 

256 In an early pioneering effort, a 3-D model derived from a terrestrial GCM was 
25? developed at the Laboratoire de Meteorologie Dynaraique (LMD). With a top at 

258 ~ 0.3 mbar, the optically active gases and haze production rates were held fixed, 

259 independent of latitude and time (Hourdin et al., 1995). Results were presented 
jeo for northern winter solstice (L s — 270°) and northern spring equinox (L s — 360°). 
26i Although the height range covered by the model limits the comparisons that can be 
362 made with the CIRS middle atmosphere results, the predicted gross seasonal changes 
263 in the 0.3 - 10 mbar region are qualitatively similar to those observed. However, the 
284 model underestimated the latitudinal temperature gradients and the strength of the 
265 winter hemisphere jet. 

ass Subsequently, Tokano et al. (1999) constructed a Titan GCM in which an attempt 
26? was made to include the effects of seasonally varying haze opacity within the frame- 
368 work of a simplified model. The potential importance of a time and spatially varying 

269 haze opacity was established, although resulting temperature gradients in latitude in 

270 the lower stratosphere were generally weaker than observed. The top of the model 

271 was at ~ 0.1 mbar so direct comparisons with CIRS results in the upper stratosphere 
2?2 cannot be made. Friedson et al. (2009) have also developed a fully three-dimensional 
273 Titan GCM, based upon the NCAR Community Atmosphere Model 3,0, including 
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in detailed calculations of surface- atmosphere interactions and radiative transfer, but 

275 with the assumption that the radiativelv active gases and aerosols were uniform in 

276 latitude and time. Their model produced stratospheric winds with a maximum am- 
577 plitude of ~ 12 ms" 1 , and equator- to-pole temperature gradients of ~ 1 K, much less 
27s than observed. 

279 More recently, the LMD group and their collaborators have developed and contin- 
5so ucd to refine a zonally symmetric model with interactive dynamics, radiative transfer, 

291 photochemistry, and cloud micro-physics (Lebonnois et al, 2001, 2003; Hourdin et ah, 
zK 2004; Rannou et ah, 2004). Horizontal transport by harorotropic eddies is parameter- 
233 ized in the 2-dimensional model using the results of Luz and Hourdin (2003) and Luz 

284 et al. (2003). Haze transported to the polar regions is found to have a major effect on 

285 the radiative heating and cooling rates resulting in enhanced net radiative cooling, 
m especially at high latitudes in the winter hemisphere. Consequently, equator-to-pole 
2 s 7 temperature gradients are increased with a corresponding strengthening of the wind 
2 ss field. Here we make use of a publiclv-available database of output from this model 
269 (Rannou et ah, 2005, http://www.lmd.jussieu.fr/titanDbase/index.hml). 

290 Results from the model have previously been compared with CIRS-derived thermal 
»i structure, zonal winds, and stratospheric gas distributions for L s ~ 300° by Crespin 

292 et ah (2008). We can now extend the comparison of the temperature and zonal wind 

293 field through L, ~ 360° with emphasis on the temporal changes. Zonal mean rnericl- 
m ional temperature and zonal wind cross sections from the model for L s — 309" and 

295 358° are shown in Pig. 10, along with the changes over this time interval. These results 

296 can be compared with the CIRS retrieved temperatures and winds and their changes 
29? shown in Figs. 2 and T. The model- predicted temperature and wind cross sections 

298 agree qualitatively with the measurements, with a strong temperature maximum in 

299 the upper levels at high northern latitudes and a strong mid latitude jet in the winter 
360 hemisphere. However, notable quantitative differences exist. The model temperature 
mi maxima for both time periods are somewhat warmer and occur at higher pressures 
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than found in the measurements. Crespin et al. (2008) suggest that this may be due 
sea at least, in part to placement of the modeled aerosol production source too deep in 
aw the atmosphere. In addition, the measurements indicate a significant increase in the 
3os stratopause altitude between the equator and the north polar region for both times, 
sot a feature that is not captured by the model. The measurements show a. temporal 
3o? decrease in temperature in the upper atmospheric levels at high northern latitudes 
3oe with the greatest change occurring at the temperature maximum centered near 0.01 
309 mbar. The model results indicate little change in the maximum with some increase 
330 in temperature at the equatorward edge of the warm region. Both the measurements 
an arid the model show an increase in zonal wind speed at, mid and low' latitudes, al- 

312 though the spatial extent and magnitude of the change differs considerably between 

313 the two. 

n* 4- 3. Potential Vorticity and Stability of the Polar Vortex 

3i i The Ertel potential vorticity and the potential temperature (a function of the 

3i6 entropy) are materially conserved quantities in adiabatic inviscid flow. In a flow that 
si? is also axisymmetric, the zonal mean angular momentum about the symmetry axis 
sis is also conserved. The zonal mean cross sections of these three quantities are thus 

319 of particular interest as tracers of the meridional circulation, or as indicators of the 

320 strength of non-adiabatic effects (radiation, for example) or of wave and eddy forcing. 

321 The seasonal changes of these three quantities should reveal the forcing that drives 

322 the general circulation. 

323 As explained in section 3.3, the retrieved temperatures can be used to make es- 

324 timates of the zonal wind. The potential temperature can be calculated from the 

325 temperature, the angular momentum from the wind, and finally the potential vortic- 

326 itv can be calculated from the wind and the potential temperature. These steps were 
22? carried out as described in Achterberg et al. (2008a). The resulting meridional cross 
326 sections are displayed in Figs. 11 through 13, as calculated from the temperatures 

and winds from Fig. 2 and Fig. 7, respectively. 
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330 Noise is magnified when the retrieved temperatures are differentiated. This has 

331 a strong effect, on the potential vorticity, because two spatial derivatives are taken. 

332 Moreover, it is the gradients of the field that are of interest in interpretations, such 

333 as baroclinic instability assessment. To reduce noise as much as possible when geu- 
33 * crating these conserved quantities, we smoothed spatially, typically with three passes 
335 of a sliding box average ten degrees of latitude in width. In addition, we evaluated 
j» temporal change by using all four time bins and fitting a linear rate of change to 

337 each (latitude, pressure) point in the meridional plane. This procedure uses all the 

338 data, and captures the important time dependence, since only a fraction of a season 
359 is available at present. 

343 The angular momentum per unit mass (Fig. 12) is dominated by the zonal mean 

341 wind, with only a minor contribution from the planetary rotation. Consequently, the 

342 angular momentum distribution is very similar to the wind field, with closed contours 

343 in the middle stratosphere. The potential vorticity cross sections (Fig, 13) show a 

344 monotonic increase from the south pole through approximately 65° N, on the poleward 

345 flank of the winter polar vortex, where there is a maximum in the potential vorticity, 

346 and then a descrease toward the north pole, where there is a local minimum. This 

347 local maximum in the potential vorticity should be viewed with caution; it is near 

348 the edge of the domain and is sensitive to the smoothing approach. Nevertheless, it 
345 consistently appears in all four of our time bins, and the amplitude of the maximum 
■iso shows an increase as time advances. 

35i There has been only limited investigation of stability criteria for nongeostrophic 
js2 vortex flows. It is known that if the Ertel potential vorticity changes sign within 

353 a vortex, or if the angular momentum per unit mass decreases with distance from 

354 the rotation axis along a surface of constant potential temperature, then the flow 

355 is unstable to symmetric or inertial overturning (Eliassen, 1951). The Titan cross 
35« sections displayed in Figs. 12 and 13 do not do this. In quasigeostrophic flow, a 
357 necessary, but not sufficient, condition for barotropic or baroclinic instabilities is a 


16 




368 change in sign of the meridional gradient of the potential vorticity (e.g. Pierrehumbert 

359 and Swanson. 1995). In the. cyclostrophic case a stability analysis has been carried 

360 out by Montgomery and Shapiro (1995), giving the same result, but the structure 
36! of perturbations is limited to low wavenumbers and therefore the conclusion is not 
363 general. The local maximum in potential vorticity would thus suggest that the polar 
363 vortex may be unstable, although it is possible that the polar vortex is stabilized by 
366 geometry or stratification. 

365 4-4- Implications for the Zonal Mean Meridional Circulation 

366 On Earth, the radiative relaxation time, r r . is roughly comparable in the tropo- 
36? sphere and middle atmosphere, approximately a few days (see e.g. Houghton, 1977; 
366 Fels, 1982; Andrews et al., 1987). In Titan’s atmosphere the range is much larger. 
36« Figure 14 shows order-of-magnitude estimates of the radiative cooling timescales at 
3?o equatorial (5°S) and north polar (75° N) latitudes, calculated using the Newtonian 
m cooling and raditation-to-space approximations (Flasar et ah, 1981; Andrews et ah, 
a?; 1987, section 2.5.2) with cooling from CH 4 , C 2 H 2 and C 2 H e , Gas abundances were 

373 assumed constant with altitude, using the CH 4 mole fraction from Niemann et ah 

374 (2005), and C 2 H 2 and C 2 H 6 mole fractions from table 2 of Coustenis et ah (2010). 
376 Because of the lack of other trace gases and of aerosols, the cooling times in Fig. 14 
m will be overestimated somewhat (see Tomasko et ah, 2008), but the differences be- 
j-7 tween the polar regions and the equator due to the temperature differences should be 
37s qualitatively correct. 

379 In the lower troposphere r r ~ 10 9 — 10 lu s (~ 100 terrestrial years or ~ 3 Titan 
3so years, Smith et al., 1981; Strobel et al., 2009). It decreases with altitude to ~ 10 7 s 
sax at 1 mbar in the stratosphere at equatorial and summer latitudes (Fig. 14) and 
382 decreases more slowly with height above that level. The slower decrease results from 
sea the compensation of the decrease in atmospheric mass per scale height with altitude, 
384 which acts to decrease r r , and the decreased emissivity of the principal hydrocarbon 
ass coolants with altitude, which acts to increase it. The variation is more extreme 
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3B6 at high latitudes in the winter hemisphere. The cold temperatures near 1 mbar 

387 imply relaxation times ~ 10 8 s. comparable to a Titan season, and near the high- 

388 latitude warm stratopause r r ~ 10 6 s. comparable to one Titan day, Hence, the 
winter stratopause region should have observable diurnal effects, but the paucity of 

sac longitude coverage from the limb scans (Achterberg efc ai, 2008a), which are needed 
m to sound these high altitudes, hamper attempts to see the variation. 

3M The large variation of r r in the middle atmosphere has implications for the struc- 
3 m ture of the zonal-mean meridional circulation. To see this, consider the zonally- 
394 averaged heat equation (e.g. Achterberg et al., 2008a): 


QT (&I g\ T~T fi , 

m~ w {dz^c p ) r r 


(2) 


393 where w is the vertical velocity, C p is the heat capacity at constant pressure, and 
is the radiative-equilibrium temperature. All variables are zonal averages. Eq.(2) 

397 says that the rate of temperature change is equal to the adiabatic heating and cooling 

398 associated with descent and ascent, plus diabatic heating and cooling, parameterized 

399 as a Newtonian relaxation of the temperature T to T eq over the time scale r,.. The 

400 balance neglects the meridional advection of heat. For cyclostrophic flow, the ratio 

401 of this to the vertical advection term in (2) is small, l/Ri, where Ri ~ 5 is the 
4G2 Richardson number in Titan’s middle atmosphere (Flasar and Conrath, 1990). In the 
403 absence of motion, (2) reduces to a purely radiative response, T r : 


d% _ Tr-Teg 

6t r T 


(3) 


404 Combining (2) and (3) yields 

djT - T r ) ~ foT l _g\ _ T - T r , 

dt 11 { dz Cp) T r 

403 Casting both the temporal change and diabatic heating/cooling in terms of T — T r 


18 
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permits a straightforward comparison of the two. The ratio of the time derivative 
*o? to the radiative term is Q s r r , where Q s is Titan’s (and Saturn’s) seasonal frequency: 
408 2 7r j 29.5 Earth years. With the exception of the cold region below the 1 mbar level at 
high northern latitudes in venter, Q s r r <C 1, and the radiative term is much larger. 
430 Hence the dominant balance in (4) is: 


w 




T-Tr 




4*3 The adiabatic, heating and cooling associated with vertical motion keeps the atmo- 
40 spheric temperature from relaxing to a purely radiative solution. This is most evident 
4i3 at high northern latitudes near the stratopause level (Achterberg et al., 2008a), where 
434 subsidence maintains a temperature, T ~ 200 K, that is much higher than the radia- 
4i5 tive solution, T r ~ 100 K (Hourdin et ah, 1995): T - T r ~ 100 K. At the 0.01 mbar 
4:6 level r r as 1 x 10® s (Fig. 14), and the static stability, (~ 4- ~ 1 K/km, taking 

437 dTjdz R: 0. From (5). this implies a subsidence of 10 cm/s, large by Titan standards. 
4i8 At low and southern latitudes, T - T r is much smaller, ~ 10 K, and r r «4x I0 6 s. 
4w Hence the relaxation term on the right-hand side of (5) is ~ 40 times smaller than 
420 near the north pole at low and southern latitudes. As the static stability at these 
423 altitudes does not vary that much with latitude, the upwelling leg of the meridional 
circulation is typified by vertical velocities < 1 nuns' 1 . Mass continuity requires that 

423 the area weighted average of the vertical velocities over latitude vanish, ignoring the 

424 possibility of hydrodynamic escape over much longer time scales (see. e.g., Strobel, 
42= 2009). It follows that, the area of subsidence near 0.01 mbar in the north is relatively 
426 compact, more than a factor of 10 smaller than the area of ascent at low and southern 
42? latitudes. 

426 Much lower in the stratosphere at high northern latitudes, Q s t t ~ 1, which is much 
423 larger than elsewhere in the middle atmosphere, and all the terms in (4) are required 
43D to account for the balance. Rannou et al. (2004) have emphasized the importance of 
aerosol transport to high latitudes in the winter hemisphere and the enhanced cooling 


43J 


432 


that can result. However, the neglect of aerosols in our estimate of r r may not he 

433 critical for the colder temperatures in the lower stratosphere, because of the extreme 

434 sensitivity of t t to temperature. Once r r becomes comparable to seasonal time scales, 

435 temperatures at high latitudes in the winter and early spring cannot decline much 
further, because the sun comes back too soon to warm up the polar region after the 

437 winter night. 

43e The 1-rnbar temperatures at 75 ° N and 85 ° N exhibit a curious behavior in late 
4» winter (Fig. 6). They increase between L s - 330 ° and 340 ° , and decrease between 

440 L$ — 340 ° and 360 ° . The increase could plausibly be explained by warming from 

441 the sun as the shadow of the polar night, recedes. The subsequent decrease, however, 

442 may reflect the interplay between solar heating and adiabatic heating from subsidence 

443 in the polar region. Evidence that the subsidence, inferred near the stratopause in 

444 the retrieved temperature cross sections (Fig. 2 ), actually extends down to the 1 
44$ mbar level comes from a consideration of the retrieved distribution of organic com- 

446 pounds aside from CH 4 in the middle atmosphere. The retrieved distribution of these 

447 gases at high northern latitudes shows enhancements down to the 1 mbar level and 

448 lower (Teanby et al., 2008 ; Coustenis et al., 2007 ). They are mostly formed from 

449 photolyt-ic, catalytic and electron-impact dissociation of N2 and CH 4 high in the at- 

450 mosphere, and they condense in the the lower stratosphere or tropopause region. The 
46i mean concentration of the organic molecules therefore increases with altitude in the 

452 middle atmosphere, and subsidence at high northern latitudes will produce higher 

453 concentrations along isobars (see e.g. Teanby et ai, 2009 ). If the 1 -mbar temperature 

454 is maintained by radiative heating and cooling and adiabatic warming from subsi- 

455 deuce, the observed cooling between L$ - 340 ° and 360 ° may signal a weakening in 
4=6 the subsiding leg of the meridional circulation at this level. If the weakening is rapid 
457 enough, it, could offset, the enhanced warming from the sun as it rises higher in the 
45s sky. 
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Table 1: Summary of Nadir Observations 


Encounter 

Start Time 

Duration 

{h:mm} 

Latitude Range 

No. of Spectra 

Resolution 
(° of arc) 

is 

TO 

'2004 Jut 02 03:30:21 

13:30 

88 ° S - 29° N 

4497 

2.1 

293 

TB 

2004 Dec 13 15:12:29 

8:25 

84° S - 67° N 

2405 

2.1 

299 

T3 

2005 Feb 14 09:57:53 

9:00 

79° S - 71° N 

2354 

2.1 

301 

T3 

2005 Feb 15 18:57:53 

4:20 

73° S - 77" N 

1585 

1.8 

301 

T4 

2005 Apr 01 08:05:16 

6:30 

76° S - 74° N 

2230 

1.9 

303 

T 6 

2005 Aug 22 20:53:37 

6:47 

83° S - 66 ° N 

2.126 

1.9 

308 

TB 

2005 Oct 27 01:24:00 

7:00 

74° S - 74° N 

3436 

2.9 

311 

TS 

2005 Oct 28 16:15:25 

7:49 

75° S . 73° N 

1915 

1.9 

31 1 

T9 

2005 Dec 27 14:04:00 

3.0:07 

74° S - 74° N 

4032 

2.9 

313 

T.tO 

2006 Jars 14 14:23:27 

9:13 

77° S - 76° N 

274.1 

2..1 

314 

Tl4 

2006 Mav 21 01:18:11 

2:00 

go _ 77 * N 

1 353 

1.7 

318 

T14 

2006 Mav 21 06:18:11 

2:58 

68 ° S - 42° N 

1900 

2.4 

318 

T15 

2006 Ju) 02 23:50:47 

7:54 

75° S - 75“ N 

1862 

2.3 

320 

TIT 

2006 Sep 06 21 .-56:61 

7:20 

66 ° S - 83° N 

2102 

2.4 

322 

T18 

2006 Sep 22 20:58:49 

7:00 

63° S - 87° N 

1705 

2.3 

323 

T19 

2006 Oct; OS 20:16:07 

6:14 

56° S - 88 ° N 

4520 

2.2 

323 

T21 

2006 Dec 11 16:02:17 

•1:03 

45° S - 88 ° N 

1248 

2.1 

326 

T22 

2006 Dec 27 15:04:13 

4:56 

38° S - 89° N 

3336 

2.0 

326 

T23 

2007 Jan 12 14:23:31 

2:15 

31° S - 52° N 

874 

2.1 

327 

T23 

2007 Jan 12 17:38:31 

2:00 

6 ° N - 89° N 

1912 

1.7 

327 

T23 

2007 Jan 13 22:38:31 

3:24 

88 ° S - 25° N 

1503 

1.9 

327 

T24 

2007 Jan 28 13:00:55 

2:15 

5° N - 89° N 

1840 

2.1 

327 

T24 

2007 Jan 29 21:15:55 

5:14 

89° S - 21 ° N 

3654 

1.9 

327 

T26 

2007 Mar 09 11:08:00 

1:41 

Sir S - 25° S 

1935 

1.8 

329 

T2T 

‘2007 Mar 25 09:07:27 

2:16 

88 ° 8 - 7° S 

2444 

1.8 

329 

T28 

2007 Apr U 12:58:00 

7:14 

46° S - 88° N 

2357 

2.2 

330 

T30 

2007 Mav 12 05:45:58 

.1:24 

88 ° S - 55° N 

1187 

1.8 

331 

T30 

2007 Mav 13 10:09:58 

1:19 

60= S - 85° N r 

1 1.23 

1.9 

331 

T32 

2007 Jun 14 07:46:11 

2:15 

2° S - 79= N 

1785 

1.9 

382 

T34 

2007 Jill 18 01:48:20 

7:23 

74° S - 73° N 

2273 

2.7 

334 

T35 

2007 Aug 13 21:32:34 

6:00 

70® S - 64“ N 

750 

2.1 

335 

T36 

2007 Oct 02 18:42:43 

8:4.6 

67° S - 72" N 

3943 

5.3 

336 

T37 

2007 Nov 19 14:47:25 

7:00 

71“ S - 81= N 

3848 

2.3 

338 

T38 

2007 Dec 05 14:06:50 

9:37 

64® S - 87® N 

6943 

2.4 

339 

T40 

2008 Jem 06 11:30:20 

7:00 

54° S - 88 ° N 

3878 

2.3 

340 

T41 

2008 Feb 23 12:32:07 

2:53 

4.6° S - 89° N 

2143 

2.7 

341 

T43 

2008 May 13 02:46:58 

6:30 

48° S - 89= N 

5255 

2.6 

344 

T44 

2008 May 27 10:24:32 

6:00 

88 ° S - 48° N 

2289 

2.5 

345 

T45 

2008 Jut 30 08:05:21 

4:07 

89° S - 49° N 

2903 

2.0 

347 

T4S 

2008 Dec 06 05:56:28 

3:06 

38° S - 89° N 

2733 

2.0 

351 

T49 

2008 Dec 20 17,24:32 

6:35 

89° S - 37° N 

3034 

2.0 

352 

T49 

2008 Dec 22 02:29:52 

3:30 

31° S - 89° N 

2884 

1.9 

852 

T51 

2009 Mar 26 11:00:31 

3:43 

89° S - 22° N 

2661 

2.0 

355 

T51 

2009 Mar 27 18:13:36 

5:12 

17® S - 89® N 

4110 

2.0 

355 

T52 

2009 Apr 03 10:29:34 

1 :48 

87= S - 18= N 

586 

2.1 

356 

T52 

2009 Apr 04 15:47:47 

7:37 

31° N - 89® N 

2427 

1.7 

356 

T54 

2009 May 05 08:11:47 

4:42 

20° S - 89° N 

4163 

3 . -.'.1 

357 

T55 

2009 Mav 21 07:09:49 

1:17 

24' S - 89" N 

.1274 

1.7 

35? 

T55 

2009 May 22 11:26:41 

8:00 

89° S - 25= N 

4264 

2.2 

357 

T57 

2009 Jun 23 08:32:35 

8:00 

89° S - 42° N 

5085 

2.2 

358 

T59 

2009 Jut 23 23:34:04 

3:00 

51° S - 89° N 

2484 

1.8 

359 

T62 

2009 Oct 1.1 15:15:21 

4:21 

70° S - 72° N 

1055 

1.9 

2 

T63 

2009 Dec 1.2 15:03:04 

5:00 

74® S - 73® N 

1996 

2.0 

4 




Table 2: Summary of Limb Observations 


Encounter 

Start Time 

Duration 

(h:mm) 

Latitude Range 

Mean 
Longitude 
(® West) 

Vertical 

Resolution 

(km) 

L s 

O 

T04 

2005 Apr 01 00:25:12 

3:30 

3° N - 86° N 

132 

35 

303 

TOG 

2005 Aug 22 01:10:00 

2:30 

34® S - 26 s N 

295 

35 

308 

TOG 

2005 Aug 22 13:40:00 

2:30 

80® S - 38° S 

318 

35 

308 

TOS 

2005 Oct 27 09:38:09 

3:20 

1 1 ° N - 86° N 

72 

40 

311 

T13 

2006 Apr 30 11:53:31 

5:00 

1® N - 41° N 

99 

28 

318 

TI6 

2006 JuJ 22 05:25:13 

2:15 

31® N - 76= N 

56 

33 

320 

T20 

2006 Oct 25 20:28:07 

3:00 

14® S - 51° N 

35 

33 

324 

T26 

2007 Mar 09 16:49:00 

4:00 

31“ S - 30° X 

305 

40 

329 

T32 

200? Jun 13 08:46:11 

4:00 

79® S - 15® N 

285 

40 

332 

T37 

2007 Nov 18 15:47:25 

4:00 

88= S- 21“ S 

158 

41 

338 

T39 

2007 Dec 21 13:57:55 

3:54 

24° S - 74® N 

99 

41 

339 

T42 

2008 Mar 25 19:27:48 

4:00 

56® S - 16° S 

179 

41 

342 

T45 

2008 Jul 31 06:58:11 

3:30 

1® S - 44® X 

190 

37 

347 

T54 

2009 May 05 13:54:10 

3:50 

21° S - 31° N 

195 

39 

357 

T59 

2009 Ju! 24 06:34:04 

4:00 

1° S - 59= N 

245 

40 

359 

T63 

2009 Dec I t 16:03:14 

4:00 

1 0 S - 83® N 

107 

39 

4 
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Figure Captions 

Figure 1. Latitude and time coverage of the nadir (top) and limb (bottom) temper- 
ature mapping observations. 

Figure 2. (top) Zonal mean temperatures, determined from both nadir and limb 
data, averaged over the time period from July 2004 through March 2006 (T0-T13 
flybys). Temperatures were averaged in 5° latitude bins, then smoothed with a 10° 
boxcar function applied three times. Contours are labeled in K. The retrieved tem- 
peratures are valid over the pressure range between about 0.005 and 8 mbar, except 
poleward of 60° N, where they are valid between about 0.005 and 2 mbar. {bottom) 
As the top panel, but averaged over the time period from July 2008 to December 
2009 (T45-T63 flybys), (middle) Difference in zonal mean temperature between the 
lower and upper panels. 

Figure 3. Northern hemispheric temperatures retrieved from all individual limb tem- 
perature maps with nearly complete northern hemisphere coverage. The dashed line 
indicates the stratopause, calculated as the pressure corresponding to the temperature 
maximum at each latitude. 

Figure 4. Equatorial temperatures retrieved from all individual limb temperature 
maps with nearly complete coverage between 30° S and 30° N. The dashed line indi- 
cates the stratopause, calculated as the pressure corresponding to the temperature 
maximum at each latitude. 

Figure 5. Southern hemisphere temperatures retrieved from all individual limb tem- 
perature maps with nearly complete southern hemisphere coverage. The dashed line 
indicates the stratopause, calculated as the pressure corresponding to the temperature 
maximum at each latitude. 

Figure 6. Retrieved 1 mbar temperatures, zonally averaged in 5° wide latitude 
bins from individual nadir maps, as a function of time, (top) Temperatures at 60° N 
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(blue squares). 70° N (red circles) and 80° N (black crosses). ( middle ) Temperatures 
as 30° S (blue squares), equator (black crosses) and 30° N (red circles), (bottom) 

Temperatures at 45° S (green diamonds). 60° S (blue squares), 75° S (red circles) and 
85° S (black crosses). 

Figure 7. Zonal winds calculated from the temperature averages in figure 2 using 
the gradient wind equation, assuming solid body rotation at the 10 mbar level at four 
times Titan’s rotation rate. Wind speed contours (black lines) are labeled in m s' 1 . 

(top) Zonal winds from temperatures over the time period from July 2004 through 
March 2006 (T0-T13 flybys). Temperatures were averaged in 5° latitude bins, then 
smoothed with a 10° boxcar function applied three times. Contours are labeled in K. 

(bottom) As the top panel, but for the time period from July 2008 to December 2009 
(T45-T63 flybys), (middle) Difference in zonal winds between the lower and upper 
panels. 

Figure 8. Least squares fits to the symmetry axis of the 1 mbar temperatures 
from nadir data from individual flybys. The left column show the fits in a sun-fixed 
reference frame, the right column in a star-fixed (inertial) reference frame. The top 
panels show the fitted azimuth of the pole fit as a function of time, the bottom panels 
show the amplitude of the offset from the IAU defined pole as a function of time. The 
dotted lines are least-squares fits to a linear variation with time. 

Figure 9. Temperatures at 0.4 mbar (red) and 1.0 mbar (blue) from Voyager 1 IRIS 
retrievals of Flasar and Conrath (1990) (points) and from CIRS data averaged over 
the period between June 2008 and June 2009 (lines). j 

j 

Figure 10. Results from a Laboratoire de Meteorologie Dynamique (LMD) 2-D circu- j 

lation model taken from http://wnyw.hnd.jussieu.fr/titanDbase/index.hmi (Rannou \ 

et a!., 2005). The three panels on the left show meridional cross sections of tempera- | 

ture for L s — 309° (top), L s = 358° (bottom), and the change in temperature between 
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the two time periods (middle). The three right-hand panels show similar results for 
the zonal wind. 

Figure 11. Potential temperature calculated from the temperature averages in fig- 
ure 2. The top panel is from data averaged from July 2004 through March 2006, the 
bottom panel is from data averaged from July 2008 through December 2009. The 
middle panel is the difference between the bottom and top panels. Contours are 
labeled in K. 

Figure 12. Zonal angular momentum calculated from the temperature averages in 
figure 2. The top panel is from data averaged from July 2004 through March 2006, 
the bottom panel is from data averaged from July 2008 through December 2009. The 
middle panel is the difference between the bottom and top panels. Contours are 
labeled in units of 10 6 m 2 s -1 . 

Figure 13. Ertel potential vortieity calculated from the temperature averages in 
figure 2. The top panel is from data averaged from July 2004 through March 2006, 
the bottom panel is from data averaged from July 2008 through December 2009. The 
middle panel is the difference between the bottom and top panels, Contours are 
labeled in m 2 K kg - 1 s _1 . 

Figure 14. Estimates of the radiative cooling timescales at 5° S (solid line) and 75° N 
(dashed line). 


32 







2005 (K) 





I 


Figure 2: 


34 



Pressure [mbar] 



Latitude 


Latitude Latitude 


Latitude Latitude Latitude 


Figure 3: 


35 




Pressure (mbar] 









Zonal Mean Temperature [K] Zonal Mean Temperature [K] Zonal Mean Temperature [K] 


300 310 


320 


330 



J§ 


I « 

i 


34Q 


It* 


350 


1 >:? 


I i i i i i I . i 1 i i I ■ i i ■ i i i i I < i I ' . 1 ,i .i-l-i .lJ-i . i. Lx^t j iXu U-xjdj-UX!LUj 



Date 


Figure 6: 


38 














Temperature (K) 






Pressure (mbar) Pressure (mbar) Pressure (mbar) 





i r - 


90 60 30 0 -30 -60 -90 

Lat 

LMD Wind Change (m s' 1 ) 


Lat 

LMD Temperature Change (K) 


Figure 10; 



2005 e 




Latitude 


Figure 11; 


S 

t 

i 


43 



Pressure (mbar) Pressure (mbar) Pressure (mbar) 



90 60 30 0 -30 -60 -90 


2009 M 



Latitude 


Figure 12: 


44 






90 60 30 0 -30 -60 -90 



90 60 30 0 -30 -60 -90 

Latitude 


Figure 13: 


45 






